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We derive the transport relaxation times for quasiparticle-vortex scattering processes via nuclear 
force, relevant for the damping of differential rotation of superfluids in the quantum liquid core of 
a neutron star. The proton scattering off the neutron vortices provides the dominant resistive force 
on the vortex lattice at all relevant temperatures in the phase where neutrons only are in the paired 
state. If protons are superconducting, a small fraction of hyperons and resonances in the normal 
state would be the dominant source of friction on neutron and proton vortex lattices at the core 
temperatures T > 10 r K. 
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A broad class of problems related to the nonequilib- 
rium spin dynamics of pulsars require the knowledge of 
the fraction of neutron star fluid interiors which is cou- 
pled to the crusts on short (unobservable) time scales. In 
practice, one needs to identify the fraction of the moment 
of inertia of the superfluid in the core of the star which 
is responding to spin perturbations as a rigid body. A 
large amount of matter in neutron star models based on 
the modern, moderately soft, equations of state (EoS) 
resides at and above twice the nuclear matter density. 
The knowledge of the state of the matter above these 
densities is limited, in particular with respect to nucle- 
onic pairing, strangeness content, etc.; thus one is led to 
consider several admissible variants of the composition of 
matter when treating the problem of transport. Here we 
shall adopt this type of an approach and examine phases 
with purely nuclconic constituents and those containing 
hypernuclear matter at large densities. 

Several electromagnetic channels of interaction be- 
tween the electrons and the nucleonic superfluids in the 
quantum liquid cores of neutron stars have been invoked 
as mechanisms of coupling of the superfluid to the normal 
matter |l]U] . Relativistic electrons couple to the anoma- 
lous magnetic moment of the neutron quasiparticles lo- 
calized in the cores of vortex lines jl| . At the densities of 
interest, the pairing in neutron matter is driven by the 
attractive P-wave interaction; a neutron vortex in a ro- 
tating 3 P2 superfluid acquires an intrinsic magnetization 
due to the pairing in the I = 1 state Q. The electron 
coupling to this magnetization dominates the coupling to 
the neutron magnetic moment at low temperatures 
If protons are superconducting even larger magnetization 
comes about from the entrainment effect |3|,[|; in addi- 
tion the proton superconductor in the mixed state might 
mediate the coupling between the superfluid neutron and 
electron fluids || . 

In this Rapid Communication we examine the nuclear 
interaction channels of coupling between the neutron vor- 
tex lines and the bulk quasiparticles; we show that (1) 
when protons are normal their scattering off the neutron 
vortex lines dominates electromagnetic coupling channels 



at all relevant densities and temperatures; and (2) if pro- 
tons are superconducting, the scattering of normal (non- 
superfluid) heavy baryons off vortices dominates the elec- 
tromagnetic coupling at the core temperatures T > 10 
K. 

To start with, consider the quasiclassical Landau- 
Boltzmann kinetic equation for the proton quasiparticlc 
distribution function 

i + ft r t ) 

dt dp dr dr dp J ' ' 

= [l-/(p,r,t)]S<(p,r,t)-/(p,r,t)E>(p ) r,t), (1) 

where s p is the proton quasiparticle energy and other 
variables have their usual meaning. The transport vertex 
in the quasiparticle limit is given by 

(n) 

S < (p,r,<) = ^ W(p'q',pq)5(£ p + E q -£ p > -e q >) 
x f(p',r,t)f(q',r,t)[l-f(q,r,t)}, (2) 



where W(p'q',pq) is the transition probability, p and q 
denote the proton and neutron momenta, respectively, 
is the neutron vortex density per unit area, and 
V is the dimensional normalization. The expression for 
scattering-out rate, E > (p, follows from Eq. (||) by 
interchanging the particle and hole occupations. 

Microscopically, the neutron vortex lattice is a collec- 
tion of spatially localized normal quasiparticles, which 
are a superposition of plain waves along the vortex circu- 
lation and two-dimensional bound states in the perpen- 
dicular direction. The continuum of Bardeen-Cooper- 
Schrieffer (BCS) paired state (neutron condensate) feels 
the intervening configuration space. The transition prob- 
ability for scattering of protons off these localized states 
(vortex cores) is given by 
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where the structure factor of the vortex core quasiparti- 
cles is expressed through the coherence factors (it, v) in 
the configuration space as 

S 1/2 (q,q') = J [u q (r ± )u q ,(r ± )-v q ,(r ± )v q (r ± )] 

x exp [i(q' x - q±) r±] dr±, (4) 

where || and _L denote components along and perpendic- 
ular to the vector of vortex circulation. The expression 
(||) relies on the Brueckner G-matrix approximation for 
the transport vertices (G R is respectively the retarded 
component) which involve the resummation of the per- 
turbation series in the particle-particle channel. 

Considering a situation not far from equilibrium we 
linearize the kinetic equation (Q) with respect to small 
departure of proton distribution function. The neutron 
distribution functions are assumed to be the equilib- 
rium Fermi distributions because the n-n relaxation time 
scales (in the relevant normal state) are much shorter 
than other time scales due to the large number density 
of neutrons. A variational method Q can be used in 
order to arrive at the relaxation time 
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v = J2^(p) 2 s(p) [WW], 
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(5) 
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ip(p) is the trial functon, and all the distribution func- 
tions are the equilibrium ones. For our purposes (macro- 
scopic momentum exchange between the quasiparticle 
current and the vortex lattice) a good trial function is 
tp(p) — P ■ Vl, where is the vortex velocity. In this 
case T> = v{p F )p 2 F v 2 L /'i(3, where v{pf) = m*p F /n 2 h 3 is 
the density of states of protons at the Fermi surface, pf 
is the Fermi momentum, m* is the effective mass of pro- 
tons, and P is the inverse temperature. Using the relation 
/(wi)[l-/(w 2 )] = g{ui-u)2)[S{ui)- /(w 2 )], where f(u) 
and gjui) are the Fermi and Bose distribution functions, 
Eq. (|5|) can be rearranged: 
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where the spectral functions of protons and neutrons are 

A<*)(k,u) = ^^|G fl (X + fc/2,ii:-fc/2;£ p+ + e p _)\ 2 
n K 
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where k = p p' and K = (p + p')/2 are the momentum 
transfer and the center-of-mass momentum, p + /_ = _KT± 
fe/2 and <?+/- = q ± fe/2. The angle intergation in Eq. 
(ph can be carried out using an angle averaged on-shell 
retarded G matrix, which then depends on the magnitude 
of the momentum transfer. To the leading order in u> /kv F 
one finds 
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\(G R (k,p F ))\ 2 9(2p F -k), (10) 



where angular brackets denote the angle averaging, v F is 
the Fermi velocity of protons, and 9 is the Heavyside step 
function. The retarded G-matrix can be further related 
to the in-medium n-p differential cross section 
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where p* n = m£m*/(m* + m*) is the reduced effective 
mass. 

The eigenfunctions of neutron core quasiparticles in 
the case of pairing (and in a spinor notation) are 
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where r, 9, z are cylindirical coordinates with the axis 
of symmetry along the vortex circulation, and p is the 
azimuthal quantum number, which assumes half- integer 
positive values. The radial functions are 
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where q± = y/q 2 — q'p, q F being the neutron Fermi mo- 
mentum, and 
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To arrive at the second equality in Eq. ( |i4[ ) we assumed 
a radial dependence of the gap function A(r) = (1 — 
e~ r /£) (Aqo being the value at r/£ — » oo), where £ is 
the coherence length of neutron condensate. In the low 
temperature limit, the transitions with a change of the 
azimuthal quantum number can be neglected, and the 
summation can be restricted to the lowest order term 
fi = fi' = 1/2. The same approximation should be used 
for the eigenvalues of quasiparticle energies. 

Substituting Eqs. in Eq. (§) and taking into 

account the normalization of states, wc find 



( k ±) = I d y Jofa v) cos ( 2x v - *■/*) e- 2K ^ 
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where Jo is the Bessel function, x = k±£ and y = r/£. 
For a fixed value of (?f£, the parametrical dependence of 
S on x can be well fitted by a Lorenzian after calculating 
the integrals numerically. For a typical value qpS, = 10, 
e.g., we find 



S(x) 



(x - x ) 2 + 7 2 



(16) 



where the amplitude is a = 0.11, the maximum point 
xq = 1.1 and the width 7 = 0.4. The neutron quasipar- 
ticle spectrum is given by 



/ih 



dr 



A (r) 2K(r) 



*U d "-""l "<{ 1+ w r h (17) 

where £° = [iirA'l /2sFn with ep n being the Fermi en- 
ergy of neutrons. The second equality is obtained by 
keeping the next-to-leading terms in small quantity q/qF 
and approximating the integral as Aqc/^. With the ap- 
proximations above the spectral function for neutrons in 
the low temperature limit (eT^ <C 1) is 



A^ n \k,oj) 



~2~ 



■exp 



-p4 



x sinh 



/2 



2 2 

w T F 

9 P 2 1,2 




(18) 



We further substitute the spectral functions in Eq. (||) 
and perform the lu integration. Further progress can be 
made if we assume that the scattering amplitude is in- 
dependent of k z and carry out the k z integration. We 
find 
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dxx S(x)—j^-(x,p F ), 
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the code of Ref. || with the separable form of the Paris 
potential and a modified Pauli-blocking operator which 
keeps only the intermediate occupation numbers of neu- 
trons and sets f(p) = (for a discussion of the related 
issues see also |j|). 

Table I shows the density and temperature dependence 
of the relaxation times r e „ (Ref. |l|), t c m (Ref. (|]), and 
T pn (present work). The microscopic parameters are from 
Ref. (l2| (the proton 1 So gap, according to |Q, vanishes 
at n = 0.4 fm , while 3 P 2 neutron gap persists through 
the whole density region of interest; recent calculations 
of Ref. which include the tensor coupling, indicate 
that superfluid neutron-normal proton mixture exists in 
the density region 0.43 < n < 0.45 fm~ 3 while matter is 
normal at larger densities). 

The ratio of electron-neutron-vortex to proton- 
neutron-vortex relaxation times T en /r p „ is of the or- 
der 10 6 and mainly reflects the difference in the spin- 
magnetic-moment and nuclear interaction matrix ele- 
ments. As the temperature decreases, both time scales 
increase and for /3 ^ £F„/7rA 2 tend to infinity due to 
the exponential cutoff of the scattering rates. Note that 
the relaxation time for electron scattering off the P-wave 
vortex magnetization, r e M, is temperature independent; 
the crossover temperature corresponding to t c m > T pn is 
of the order of 3 x 10 6 K. 

Consider now neutron star matter which, in addition 
to the neutron-proton-electron (npe) mixture, contains 
also heavy baryons in the normal state. Among these 
the lowest threshold densities have the E~ and A hy- 
perons and the A - resonance. The threshold densities 
in the case of noninteracting Fermi gases and relativistic 
mean-field models are low, a few times nuclear saturation 
density (n sat = 0.17 fm~ 3 ) fiq , |ll|| . Correlations tend to 
shift the threshold densities to higher values in general. 

TABLE I. The relaxation times for the spin period of the 
Vela pulsar P = 0.089 s (r scale oc P) and for temperatures 
~ 10 s , 10 7 , and 3 x 10 6 K (in descending order). x v is 
the proton fraction. 



with x\ = 2pf£,/?i; here Kq is the modified Bessel func- 
tion. The temperature dependence of the relaxation time 
is due to protons contributing the (thermal smear- 
ing of the Fermi surface) and vortex core quasiparticles 
contributing the exp(— /3e^ 2 ) (the probability of excit- 
ing a neutron quasiparticle inside the core of a neutron 
vortex). The n-p scattering cross section is modified in 
the medium, first, due to the suppression of the inter- 
mediate state two-particle propagation and, second, due 
to the modification of the density of states. We incorpo- 
rate the medium modifications of the cross section using 



n 

(fm- 3 ) 


.i p 


A 

(MeV) 


[fm 2 ] 








0.4 


0.060 


0.7 


87.2 


0.001 s 
0.004 s 
3.670 s 


0.027 d 
0.108 d 
11.190 d 


87.175 d 


0.6 


0.070 


1.3 


51.7 


0.006 s 
0.015 s 
41.56 d 


4.516 d 
12.426 d 
18947. 12 yr 


51.699 d 


0.8 


0.084 


1.3 


57.7 


0.011 s 
0.022 s 
111.71 d 


9.791 d 
19.175 d 
36279.45 yr 


57.865 d 


1.0 


0.090 


1.2 


65.1 


0.015 s 
0.024 s 
23.99 d 


6.876 d 
10.811 d 
397320 yr 


65.098 d 
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While the uncertainties in the hypernuclear matter cal- 
culations are still considerable, we shall, for the purpose 
of order of magnitude estimate, make several simplifying 
assumptions below. In the high density n^peY phase, 
where superscript s stands for the superfluid state and 
Y = £~, A, A - , the Y-n channel due to the heavy 
baryons would contribute additively to the scattering 
rate in the p-n channel; in this case no new qualitative 
effects are expected. 

The situation is different in the n^p^eY phase, 
where in addition to the neutron condensate, the pro- 
ton superconductor is in the mixed state. The relaxation 
time for the dominant electromagnetic coupling via scat- 
tering of electrons by the magnetic flux of a vortex is 
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where k e p is the electron Fermi wave vector and A is the 
penetration depth of the proton superconductor, ru> 
is the density of either neutron or proton vortex lines. 
To compare this time with the case of nuclear scat- 
tering of heavy baryons by neutron and proton vortex 
core quasiparticles we assume in Eq. (|l^) (with obvious 
modifications) oy„ = 40 mb and replace the effective 
masses of heavy baryons by their bare masses mjy = 1116, 
m s - = 1197, and m A - = 1232 MeV. Once the threshold 
density is achieved, the hyperon population grows steeply 
to about 10% and stays further almost constant Jl0| , pl| . 
Table II shows the net relaxation time [Eq. (|n])] due to 
A , £ _ , A - scattering off neutron or proton quasiparti- 
cles localized in the vortex cores at the threshold density 
n = 0.4 fm -3 , and for a well developed hyperon core 
at n = 0.5 fm -3 . The relative abundances for heavy 
baryons are assumed equal and fixed at 1% and 10% for 
each case respectively. The values of Ty v are larger than 
that of TYn at n = 0.4 fm~ 3 mainly due to the large ratio 
of neutron to proton Fermi energies, i.e., more abundant 
neutrons are more effective in scattering of hyperons. 

TABLE II. Relaxation times in hypernuclear matter; con- 
ventions are the same as in Table I. xy is the concentration of 
each species of heavy baryons (£~, A, and A~); TYn and ty p 
are the relaxation times for scattering of all species of heavy 
baryon off neutron and proton vortex core quasiparticles, re- 
spectively. The electron-vortex-flux relaxation time, r e $, is 
given per single flux quantum <f>o — nhe/c. 



n 

(fm" 3 ) 




x Y 


A p 
(MeV) 


A n 
(MeV) 


TYp 


TYn 














0.036 s 


0.0002 s 




0.4 


0.05 


0.01 


0.4 


0.7 


70.34 s 


0.024 s 


0.167 s 












292.4 d 


21.18 s 




0.5 


0.1 


0.1 


0.1 


1.0 


0.027 s 
0.32 s 
1.43 s 


0.0002 s 
15.15 s 
45.94 d 


0.195 s 



At n = 0.5 fm~ 3 the trend is opposite mainly due to an 
order of magnitude difference in the pairing gaps. One 
should note, however, that for typical pulsar magnetic 
fields, B ~ 10 12 — 10 13 G, the number of proton vortices 
per neutron vortex is very large, /n^ ~ 10 12 -10 13 
and Typ is the main quantity of interest. For core temper- 
atures T ~ 10 8 K, typical for not extremely old neutron 
stars, the hyperon scattering appears to be the dominant 
source of friction. It could be also important in the well 
developed hyperon cores at T ~ 10 7 K, below the den- 
sities at which the proton superconductivity is quenched 
because of the high proton fraction. 

In summary, our calculations above show that scatter- 
ing of baryons via nuclear force off vortex core quasi- 
particles in neutron star's superfluid interiors may be 
dominant in a number of cases. In the so-called nS s 'pe 
(or n^peY) phase, where only the neutrons form a con- 
densate, the baryon (including hyperons and resonances) 
scattering off neutron vortex core quasiparticles domi- 
nates other channels driven by the electromagnetic inter- 
actions [0J|]. In the n^p^e phase, where the type-II 
proton superconductor is in the mixed state, a trace of 
heavy baryons (such as £ _ ,A, and A - ) in the normal 
state would produce a resistive force on vortices which 
will dominate at temperatures ~ 10 s K and will give 
the way to electron flux scattering at lower temperatures 
< 10 7 K Which of the channels dominates at the 

intermediate temperatures ~ 10 7 K, depends on the de- 
tails of the microphysical input, such as the pairing gaps 
and baryon abundances. For practical applications, our 
main result, Eq. ([i"9|), can be cast in a compact form 
(fi = c= 1) 



Tpn — f -5 



pn 



kF P 

kFn 



13 expf 0.038 )P, (22) 



pn 



ffi>n k% 



where the relations = (4m n /?i)p- 1 , where P is the 
spin period, and ct*„/47t = da* n /dQ, have been used. 

As an illustrative example consider first a neutron star 
model of canonical mass 1.4 M© based on the EoS of Ref. 
|l3| . The n^pe-phase in this model occupies the radial 
range < R < 8 km for the gap profiles of Ref. |l2| and 
7 < R < 7.7 km for the profiles of Ref. |l4| (star radius is 
10.13 km). For the maximum mass, M = 1.65Af Q , star 
model of Ref. |ll| (star radius is 11.4 km), the n^p^eY- 
phase occupies the shell 7 < R < 9.5 while at smaller 
radii the 1 So proton pairing vanishes (n^peF-phase). 
These numbers by no means provide a consistent overall 
picture of the distribution of different superfluid phases 
in a neutron star, for the results on the pairing gaps 
and EoS stem from different (BCS-BHFj||Jl|l vs respec- 
tively variational 1 13 1 and mean-field |10[|) approaches to 
the nuclear many-body problem. Nevertheless, our obser- 
vation of the dominance of nuclear interaction channels 
when any kind of baryons is in the normal state is ro- 
bust, i.e., independent of the details of the microphysical 
input. Its implementation, however, as given by the ex- 
ample above, is connected with the uncertainties due to 
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limited knowledge of the interactions and pairing in the 
high density nuclear matter as well as the heavy baryon 
concentrations and their density thresholds. 

This research has been supported by the Max Kade 
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